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Abstract 

In this paper propagation properties of a parallel-plate waveguide with tunable artificial impedance surfaces as 
sidewalls are studied both analytically and numerically. The impedance surfaces comprise an array of patches over 
\^ • a dielectric slab with embedded metallic vias. The tunability of surfaces is achieved with varactors. Simple design 

equations for tunable artificial impedance surfaces as well as dispersion equations for the TE and TM modes are 
' ^ ' ■ presented. The propagation properties are studied in three different regimes: a multi-mode waveguide, a single-mode 

Q^l waveguide, and below-cutoff waveguide. The analytical results are verified with numerical simulations. 
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I. Introduction 

Artificial impedance surfaces [l]-[9] have received a lot of interest since the beginning of the last 
decade. In general the artificial impedance surfaces are composed of a capacitive grid on top of a thin 
metal-backed dielectric substrate. The substrate may include vias, as in [2], [4], or may not include them 
psj ( [9]-[ll]). Nevertheless, the purpose in both types of designs is to use the grounded substrate to provide 
> an inductive response that together with the capacitive grid would create a resonant structure. Because of 
^ the resonant nature such impedance surfaces are commonly referred to as high-impedance surfaces. 
^ Recently some research has been devoted to electrical tunability of high-impedance surfaces. The 
tunability has been realized e.g. by connecting adjacent patches or strips to each other by voltage 
0\ • controllable varactors in order to construct tunable antennas [12], [13], phase shifters [14], lenses [15], 
^ band-pass filters [16], and band-stop filters [17]. In addition to the varactor-based tunable impedance 
\ surfaces also a MEMS -based tunable high-impedance surface has been proposed in [18]. In this paper we 
^ \ concentrate on the varactor-based tunable high-impedance surfaces. 

In [12]-[15] the analysis of the tunable high-impedance surface has been done by using a simplistic 
lumped-element model of the surface, similar to that of [19]. In addition to the lumped element models, 
5^ also a layered homogeneous material models have been developed [20], [21] to predict the behavior of 
non-tunable high-impedance surfaces. The lumped-element model [19] offers guidelines for the design of 
a tunable high-impedance surface. For instance, in [22] an equivalent circuit model is used to approximate 
the effect of the varactor resistance to the reflection characteristic of a high-impedance surface similar to 
that described in [12]. However, the lumped-element model cannot be used for an accurate analysis of 
the surface nor for the analysis of many applications because it does not take the oblique incidence into 
account. 
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Because of their unique characteristics the artificial impedance surfaces have been used as wall coatings 
in different waveguiding structures. For instance, close to the resonance frequency the input surface 
impedance of an artificial impedance surface is high and the surface behaves as a magnetic conductor. 
This feature has been utilized in quasi-TEM waveguides [11], [23]-[25]. Furthermore, the possibility to 
electrically vary the input surface impedance of an artificial impedance has been exploited in many waveg- 
uiding applications [14]-[17]. For the design of such applications accurate knowledge about the waveguide 
modes is essential. In [26] analysis of dispersion in a rectangular waveguide with impedance sidewalls 
comprising non-tunable dipole-like frequency selective surfaces (FSS) on a metal-backed dielectric slab 
has been done numerically by using the method of moments. 

In this paper we introduce a simple analytical model for a varactor-tunable high-impedance surfaces that 
predicts the response of the impedance surface very well even for oblique incidences. This model is general 
and can be used for surfaces that comprise any type of a rectangular patch array, for instance for such as 
those in [12], [13]. Together with the dispersion equations, the model for the tunable impedance surfaces 
is used to study the propagation properties of a parallel-plate waveguide having either one or two tunable 
impedance surfaces replacing the metallic plates. In particular, parallel-plate waveguides allowing multi- 
mode or single-mode propagation are studied. Also a waveguide operating below its cutoff frequency is 
considered. Numerical full-wave simulations verify the analytical results and show that the used analytical 
model describes accurately the properties of tunable impedance surfaces in waveguiding set-ups. 



We study the propagation properties of a parallel-plate waveguide having tunable impedance surfaces 
using the plane-wave interpretation. The two-dimensional waveguide geometry is illustrated in Fig.[T] The 
number of the transversal wave numbers is limited to one as no propagation takes place in the x-direction. 
The waveguide is confined in the ^/-direction by plates that can be modeled with an impedance surface Z^^ 
or Zj^p that are dependent both on the wave number k and the propagation constant along the waveguide 
p. The notation ± refers to the upper/lower surface, respectively. 

The dispersion relations for the parallel-plate waveguide presented in Fig. \T\ can be solved from the 
boundary conditions at the waveguide sidewalls and from the transverse components of the magnetic and 
electric fields. The transverse components can be calculated from the longitudinal (z-) components using 
general plane-wave solutions in the y-dimension (see e.g. [27]). The resulting dispersion equation for the 
TE modes can be written as: 



II. Analytical model 



taii{kyd) = jr] 




(1) 



For TM modes the dispersion equation reads 
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Fig. L Illustration of the two dimensional waveguide confined by two impedance surfaces. 
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(a) (b) 
Fig. 2. (a) A view from the top of the Sievenpiper mushroom structure loaded with varactors. (b) A view from the side. 



Fig. 3. A transmission-line model for modeling artificial impedance surfaces. Zq is the free-space impedance. 

In the above formulae ky = ±^/k'^ — /S"^ is the transverse wave number, and r] is the plane-wave impedance 
of the medium filling the waveguide. By choosing the "minus"-branch of the transverse wave number, 
also surface modes on the waveguide sidewalls are predicted by the above dispersion equations. 



A. The input surface impedance 

The dispersion equations © and Q have been derived for arbitrary surface impedances. In this paper the 
tunable version of a mushroom-type artificial impedance surface (see Fig.[2l), proposed in [2], is studied as 
a possible particular realization. The metallic plates of a waveguide are replaced with artificial impedance 
surfaces that are comprised of metallic rectangular patches and metal-backed dielectric substrates with 
embedded vias. The tunability is achieved by connecting the adjacent rectangular patches to each other 
by varactors. The input surface impedance of an artificial impedance surface can be modeled through a 
transmission-line model shown in Fig. [3l The input impedance is hence a parallel connection of the grid 
impedance of an array of rectangular patches and the surface impedance of a metal-backed dielectric slab 
with embedded vias: 



inp 



(3) 



In the above equation subscript g refers to the grid impedance of an array of rectangular patches and s 
refers to the surface impedance of the substrate. 

A simple and accurate analytical model for a non-tunable mushroom-type impedance surface is available 
[7]. In [7] the mushroom structures comprised arrays of patches on top of a dielectric layer. However, 
in this paper we consider mushroom structures that comprise arrays of patches on top of a metal-backed 
dielectric slab with embedded metallic vias. The vias are needed to provide the bias voltage for the 



4 



varactors that are used to vary the capacitance between the adjacent strips or patches (as in [12], [14]). 
For TE modes the electric field component is perpendicular to the vias and, in the case of thin vias, does 
not excite them. It can be concluded that the analytical model [7] can be readily applied for TE modes 
in the case of embedded vias by taking the effect of varactors into account correctly. However, for TM 
modes the electric field has a parallel component to the vias. In this case it is possible to model the 
metallic vias in a dielectric slab as an effective wire medium [3], [5]. 



B. Surface impedance of the dielectric slab with embedded vias 

The surface impedance for a wire medium comprising thin perfectly conducting wires reads [5]: 



jw/io 7^ 7T-^, (4) 

7tm - 



where k = ko^/e^ is the wave number in the host medium. 



7tm = i^'^SQEtHQ - — /3^ (5) 



Ex. is the relative permittivity for the fields along the transverse plane, and 

1 



K = , (6) 
J- In 

2-iT 4ro(a— ro) 

Furthermore, a is the period of the wires, ro is the radius of the wires, and the relative permittivity for 
the fields along the normal of the medium reads 

e^ = eAl- ^\ . (7) 



In the case when the vias are thin and vertically oriented, the relative permittivity for the fields along the 
transversal plane, St, equals to the relative permittivity of the host medium, Sr. 

For TE modes the electric field is perpendicular to the thin metallic wires. In this case, as discussed 
above, the electric field does not excite the wires, and the surface impedance for the TE mode is that of 
a metal-backed dielectric slab (see, e.g., [5]): 

„TE tan{kyy/e^h) 

— J^l^o — 7 — 1= — • yo) 



C. The grid impedance 

The grid impedance for an array of patches can be calculated through the approximative Babinet 
principle using the averaged boundary conditions for a mesh of wires or strips. The averaged boundary 
condition for a mesh of strips is available e.g., in [5]. The grid impedance for an array of ideally conducting 
patches on top of a dielectric substrate reads [7]: 

z™ = (9) 



where the effective wave impedance rics = -7^, the effective wave number A^cfr = kQ^Je^, sin (9) = 

V^cff * 

^ , and b and a are the dimensions of the unit cell of the structure along x— and 5;— axis, respectively. 
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Further, a is tlie grid parameter. 



- = ^ln|-^l, (11) 
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sin(l^) 



where w is the gap between the adjacent patches (see Fig. [2l). The limitations for above expressions of 
the grid impedance have been discussed in more detail in [7]. It can be concluded here that ^ and (flOl) 
are valid when w < a,b and up to the frequencies when a,b ^ |. The effective relative permittivity for 
the array of patches or grid of strips on the boundary between two medium having relative permittivities 
of El and 62 reads approximately [28]: 

Scs = — ^ — • (12) 

In this paper the array of patches is located on a boundary between free space and wire medium. From 
above it is known that in a wire medium the fields along the transversal and normal axis see different 
effective relative permittivities. For an array of patches the electric fields are concentrated mainly between 
the adjacent patches, transverse with respect to the vias, and the effect of the vias on the electric response 
is weak. For this reason the transversal relative permittivity of the wire medium is used in (fT2l) . 

The grid impedances Q and (flOl) can be written in a lumped-element form as: 

^TM,TE _ 1- /I ox 

g ~ • ^TM,TE' ^ 

where C™''^^ is the grid capacitance for the TM- or TE-polarized incidence fields. The grid capacitance 
equals to the averaged capacitance per one unit cell in the x— and direction for the TM- and TE- 
polarized cases. Using ([9]), (flOl) . and (fT3l) the grid capacitance for an array of patches can be written 
as: 

l^sin(^)'' 

^r-^^^^^mf^Vi-f^-l- (15) 



TT 



The above formulas for the capacitive grid impedances hold for ideally conducting patches. Although 
the capacitive impedance is derived through averaged fields on the grid, we may consider a lumped- 
element model here. This way the capacitance of the varactors, that are used for tuning, can be included 
in the analysis easily: The additional capacitance of the varactors is connected in parallel with the grid 
capacitance and the total impedance of a unit cell is thus a parallel connection of these impedances. 
Hence, the total grid capacitance of the array of patches with varactors can be written as: 

^TM,TE _ ^TM,TE . r< ^^f.\ 
^tot ~ '-^g "-"vari 

where Cg is the grid capacitance for an array of ideally conducting patches and Cvar is the capacitance 
of the tunable varactor. Finally, the total grid impedance reads: 

7TM,TE _ 1 _ \ /i-yN 

^g ~ . ^TM,TE ~ / \ • V^'-' 



In order to most effectively tune the grid impedance, the capacitance of the varactor needs obviously 
to be considerably larger than the grid capacitance. 

The effect of the varactor resistance to the performance of a tunable high-impedance surface has been 
studied in [22]. It is possible to include the effect of the intrinsic diode resistance to the above analysis 
in a similar way as done in [22]. However, this is considered to be out of the scope of this paper. 



6 




Fig. 4. Color online. Simulation model of the tunable high-impedance surface unit cell. Because of the periodicity, the capacitance of each 
capacitive sheet equals to 2Cvar- 
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Fig. 5. Color online. The reflection phases for different incidences and for different values of varactor capacitance. HFSS simulation results 
are denoted with crosses and circles for varactor capacitance values of 60 fF and 120 fF, respectively. The simulation results for normal 
incidence, 45° (TM), and 45° (TE) are colored blue, green, and red, respectively. 

III. Numerical validation and results 

In this section the propagation properties of a parallel-plate waveguide having either one or two artificial 
impedance surfaces are studied analytically. The analytical results are verified with simulations using 
Ansoft's High Frequency Structure Simulator (HFSS). 

The impedance surfaces for the waveguides are designed for the lower millimeter wave region, namely 
Ka-band (26-40 GHz). Following the notations in Fig.[2l the parameters of the studied artificial impedance 
surface are: a = b = 1 mm, w = 0.1 mm, h = 0.2 mm, and £r = £2 = 4. The medium inside the impedance 
waveguide is air (ei = 1). As in any resonant circuit, the bandwidth of the resonance can be increased by 
increasing the effective inductance (the height of the grounded dielectric slab) with respect to the effective 
capacitance. Here the resonance bandwidth of the high-impedance surface is decreased intentionally so 
that the effects due to the resonance are clearly recognizable in the dispersion plots. The varactors on 
each edge of the patch are modeled with lumped capacitive sheets whose value of capacitance is changed 
depending on the studied case. The simulation model of the high-impedance surface is shown in Fig. HI 
The periodicity in the simulation model is achieved by using the periodical boundary conditions available 
in HFSS. 

According to (fT4l) the grid capacitance of the designed impedance surface is approximately 26 fF. Based 
on this information and knowing the frequency band of interest, the capacitance of the varactors is varied 
from 60 fF to 120fF; with these values the resonance frequency of the surface for the normal incidence 
appears to be approximately at 34 GHz and 26 GHz, respectively. Furthermore, the capacitance of the 
varactors is considerably larger than the grid capacitance, as discussed earlier. 
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Fig. 6. Color online, (a) The propagation properties for TE modes in an impedance waveguide with one tunable impedance surface, (b) The 
propagation properties for TM modes in an impedance waveguide with one tunable impedance surface. The fundamental modes of metal 
waveguide are plotted with dash-dot lines. /3te and /3tm refer here to the propagation constants of the TE and TM modes, respectively. 
Only the two lowest modes are simulated. 



The reflection phases of the surface for different values of Cvar and for different incident angles are 
shown in Fig. [5l The analytical results concur with the simulation results very accurately. Clearly the 
bandwidth for the TE polarization becomes smaller and the bandwidth for the TM polarization larger as 
the angle of incidence grows. The effect of this to the propagation properties of the impedance waveguide 
will be discussed later. 

A. A multi-mode waveguide 

Before moving on to the interpretation of the dispersion plots, the terms used in this and following 
paragraphs need to be clarified. It is natural to use the modes propagating in an empty metallic parallel- 
plate waveguide as reference cases when studying the propagation properties of an impedance waveguide. 
Therefore the modes propagating in a metallic waveguide will be referred as fundamental modes from 
here on in order to distinguish them from the modes of an impedance waveguide (referred plainly as 
modes). 
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(a) (b) (c) 

Fig. 7. The normalized magnitude of tlie a-— component of the electric field for the TEoi mode at points (a) fib = 0.5, / = 25.0 GHz (b) 
/3b = 0.69, / = 33.5 GHz, and (c) /36 = 0.71, / = 34.0 GHz 



The dispersion curves of a 20-mm high parallel-plate waveguide with one tunable impedance surface 
are shown in Fig. [6l The second surface is perfectly conducting metal. The height of the waveguide was 
chosen so that many fundamental modes would propagate in the waveguide at Ka-band. The fundamental 
modes of a parallel-plate waveguide are shown with dash-dot lines. The simulation results have been 
marked by crosses in Fig. [6l The concurrence between the analytical and numerical results is very good. 

In Fig. |6] near the resonance frequency of the impedance surface a 180° mode hop is observed for 
both polarizations. The propagating wave exhibits a 180° phase shift while the mode morphs from one 
fundamental mode to another. This is due to the change of the reflection phase of the impedance surface. 
For instance, a wave propagating at 27.5 GHz near the second fundamental mode, TE02 in Fig. [6ta), 
morphs into the first fundamental mode, TEqi, when the varactor capacitance is changed gradually from 
60 fF to 120fF over a certain distance. This is similar to the behavior of corrugated-waveguide mode 
converters, where the depth of the corrugations is tapered gradually [29]. Similar behavior occurs for 
higher-order modes and for TMqi as well. The field pattern of the TEqi mode is plotted in Fig. U\ 

The mode conversion is seen to happen more gradually for the TM modes in Fig. [6](b) than for the TE 
modes in Fig. [6] (a). This is because the bandwidth of the high-impedance surface becomes wider for TM 
polarized fields than for TE polarized fields as the angle of incidence grows, as discussed above. This 
creates an advantage for the TE modes over the TM modes in tunable impedance waveguide applications: 
The needed range of tuning of the resonant frequency of the high-impedance surface is smaller for the 
TE modes than for the TM modes in impedance waveguide applications. 

The dispersion curve for a parallel-plate waveguide having two tunable impedance surfaces is shown 
in Fig. [H Similar mode conversion is seen as in the case of just one impedance surface. However, instead 
of a 180° mode hop discussed earlier, a 360° mode hop occurs. This is simply because both impedance 
surfaces induce a 180° mode hop. 

In the case of two tunable impedance surfaces (see Fig. [8] (a)) interesting properties in the vicinity of 
the resonance frequency of the surface are seen: Both TEqi and TE02 modes cross the light line at only 
slightly different points. The TEqi mode crosses the light line at (3b = 0.71, / = 33.9 GHz whereas the 
TE02 mode crosses it at f3b = 0.77, / = 34.2 GHz. The field patterns of the TEqi and TE02 modes are 
shown in Figs. |9]and[10l respectively, for the case when Cvar = 60 pF. Clearly, the modes can be divided 
into symmetric and asymmetric modes which is not possible in the case of just one tunable impedance 
wall. 

B. A single-mode waveguide 

The dispersion curves of a 7-mm high waveguide having one tunable impedance- surface are shown 
in Fig. \TT\ (a) and in Fig. \TT\ (b) for a waveguide having two tunable impedance- surfaces. In a regular 
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Fig. 8. Color online, (a) The propagation properties for TE modes in an impedance waveguide with two tunable impedance surfaces, (b) The 
propagation properties for TM modes in an impedance waveguide with two tunable impedance surfaces. The fundamental modes of metal 
waveguide are plotted with dash-dot lines. /3te and /3tm refer here to the propagation constants of the TE and TM modes, respectively. 
Only the two lowest modes are simulated 



metallic waveguide only one mode would propagate in the waveguide in the Ka-band. However, for the 
impedance waveguide, Fig. [H] clearly shows two (three) modes in the case of one (two) high-impedance 
sidewall(s). Because the response of the high-impedance surface for the normal incidence is the same for 
both TE and TM modes, also the cutoff frequencies are the same. 

In Fig.[ll](a) the first TM mode diverges from the first fundamental mode and converges to the light line 
as we move through the resonance frequency of the high-impedance surface up to higher frequencies. The 
first TE mode diverges also from the first fundamental mode and crosses the light line at / ^ 34 GHz (the 
resonance frequency of the high-impedance surface having C^ar = 60 fF) after which the mode transforms 
into a surface-wave mode. The second mode in the waveguide has the cut-off frequency of / = 33.3 GHz. 
Both of these modes converge gradually to the first fundamental mode. 

In the case of two impedance surfaces (Fig. \TT\ (b)) we can see two symmetric modes with cut-off 
frequencies of 19.7GHz and 35.3 GHz and one asymmetric mode with a cut-off frequency of 31.6 GHz. 
Both the first symmetric and asymmetric TE modes cross the light line at / = 33.8 GHz and 34.8 GHz, 
respectively, and transform into surface-wave modes. The first symmetric TM mode crosses the light line 
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Fig. 9. The normalized magnitude of tlie a-— component of the electric field for the TEoi mode at points (a) pb = 0.5, / = 24.9 GHz (b) 
/3b = 0.69, / = 33.3 GHz, and (c) /36 = 0.71, / = 33.9 GHz 
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Fig. 10. The normalized magnitude of the x— component of the electric field for the TE02 mode at points (a) j3b — 0.5, / — 27.7 GHz 
(b) Pb = 0.69, / = 33.9 GHz, and (c) f3b = 0.77, / = 34.2 GHz 

at / = 27.8 GHz and the second symmetric TM mode converges to the first fundamental mode. The 
asymmetric TM mode converges to the light line. This means that there exists a stop band for both TM 
and TE modes. For TE modes the width of the stop band is 1.5 GHz between the symmetric TE waveguide 
modes and for the symmetric TM modes the width equals 7.5 GHz. 

C. A below-cutoff waveguide 

In Figs.[T2](a) and (b) the dispersion curves of a 3.5mm-high parallel-plate waveguide are shown for the 
cases of one and two impedance sidewalls, respectively. Fig. [12] shows that although no TE mode would 
propagate in the metallic waveguide below 42.8 GHz, one mode propagates in the impedance waveguide 
with one impedance sidewalls, and two modes propagate in the case of two impedance sidewalls. Similarly, 
in the case of TM modes, we find one or two extra modes below the cut-off frequency depending on the 
number of impedance sidewalls. 

The TE modes form a narrow pass band both in the case of one and two impedance sidewalls. In the 
case of just one impedance sidewall the cutoff frequency of the first mode is 31.6 GHz and the TE mode 
crosses the light line at 34.8 GHz. In the case of two impedance sidewalls the cutoff frequencies for the 
symmetric and the asymmetric modes equal 29.5 GHz and 34.5 GHz, respectively. The symmetric mode 
crosses the light line at 33.9 GHZ and the asymmetric mode crosses the light line at 35.8 GHz. 

In the case of two impedance sidewalls both propagating TM modes demonstrate backward-wave 
propagation (see Fig. [T3] ). The symmetric TM mode (cutoff frequency at 29.5 GHz) first propagates as 
a forward wave but transforms into a backward wave after the point Pd = 0.4. The asymmetric mode 
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Fig. 11. Color online. The propagation properties of an 7-mm high impedance waveguide (a) with one tunable impedance surface and (b) 
with two tunable impedance surfaces. The value of the varactor capacitance is Cvar = 60 fF. The fundamental modes of metal waveguide 
are plotted with dash-dot lines. /3te and /3tm refer here to the propagation constants of the TE and TM modes, respectively. 

(cutoff frequency at 34.5 GHz) is first a backward wave and transforms to a forward wave after the point 
Pd = 0.35. In a 3.5 -mm metallic waveguide a TEM mode having the orientation of the field components 
similar to the considered TM mode would propagate. Figs. [Ill (a) and (b) show that both with one or two 
tunable impedance surfaces it is possible to create a tunable stop band for the TEM mode. 

IV. Discussion and conclusions 

An analytical model for a general type of tunable impedance surface comprising an array of patches has 
been presented. Dispersion equations for a parallel-plate waveguide having arbitrary surface impedance 
sidewalls have been presented. Together with the presented dispersion equations, the analytical model for 
the tunable high-impedance surfaces is used to study the propagation properties of impedance waveguide 
having either one or two tunable impedance surfaces. In oversized waveguide mode conversion is shown. 
In single-mode waveguide multi-mode propagation and band gaps are shown. Furthermore, forward- as 
well as backward-wave propagation in a below-cutoff waveguide is presented. The results are validated 
with commercial full-wave simulation software. The concurrence between the analytical and numerical 
results is very good. 
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Fig. 12. Color online, (a) The propagation properties of a below-cutoff impedance waveguide with one tunable impedance surface, (b) The 
propagation properties of a below-cutoff impedance waveguide with two tunable impedance surfaces. The value of the varactor capacitance 
is Cvar = 60 fF. /9te and /Jtm refer here to the propagation constants of the TE and TM modes, respectively. 



The presented analytical model for the study of the propagation properties in an impedance waveguide 
has proven to be very useful. The model has been verified and used to predict the dispersion in tunable 
impedance-wall waveguides in various example cases. The time needed to produce the results with the 
analytical model compared to the time consumed by numerical simulations is very marginal. 

The dispersion results of the impedance waveguide show interesting features. In oversized waveguides 
tunable impedance surfaces allow one to realize tunable mode converters [29] and in many other applica- 
tion, where field transforming inside waveguides is needed. In addition, in single-mode or below the cutoff- 
frequency waveguides possible applications for the tunable impedance surfaces include phase shifters [14], 
filters [16], [17], and lenses [15]. Without doubt, the design work for all these applications benefit from an 
efficient and simple model to predict the dispersion characteristics of an impedance waveguide. This list 
can be continued with such applications, as different types of antennas [12], [13], [30], tunable artificial 
magnetic conductors (AMC) [9], [31], and tunable electromagnetic band-gap structures (EBG) [9], [32], 
[33]. 
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Fig. 13. A zoom-in of the Fig. [12] The backward wave propagation of the both symmetric and asymmetric TM modes. 
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